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Boundary Effects on Dynamic Behavior of Josephson-Junction 

Arrays 



The boundary effects on the current-voltage characteristics in two- 
dimensional arrays of resistively shunted Josephson junctions are examined. 
In particular, we consider both the conventional boundary conditions (CBC) 
and the fluctuating twist boundary conditions (FTBC), and make comparison 
of the obtained results. It is observed that the CBC, which have been widely 
adopted in existing simulations, may give a problem in scaling, arising from 
rather large boundary effects; the FTBC in general turn out to be effective 
in reducing the finite-size effects, yielding results with good scaling behavior. 
To resolve the discrepancy between the two boundary conditions, we propose 
that the proper scaling in the CBC should be performed with the boundary 
data discarded: This is shown to give results which indeed scale well and are 
the same as those from the FTBC. 
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During the past decades superconducting arrays have been extensively studied both 
theoretically and experimentally [IJ. The superconducting arrays in equilibrium are well 
described by the XY model, which exhibits the Berezinskii-Kosterlitz-Thouless (BKT) tran- 
sition at a nonzero temperature in two dimensions 0. Theoretically this transition is man- 
ifested by the characteristic behavior of static quantities such as the universal jump of the 
helicity modulus. On the other hand, the standard experimental probe of the transition 
is to measure the dynamical response to the external current driving |$|. The need for 
the analysis of the experimental data together with the interest in the dynamic properties 
themselves has motivated extensive theoretical examinations of the dynamics of the system, 
where numerical simulations have been of great help. 

It has been recognized that the boundary conditions employed in those simulations, per- 
formed in finite systems, may affect the results in a crucial way, which makes it necessary 
to examine their effects. Most existing dynamic simulations of the system in the absence of 
driving currents have been performed with periodic boundary conditions along both direc- 
tions; under driving currents, the boundary conditions along the current-injected direction 
are changed to be free whereas periodic boundary conditions are retained along the per- 
pendicular direction. Such conventional boundary conditions (CBC), however, are known 
to produce rather large edge effects, and the system is expected to approach the thermo- 
dynamic limit rather slowly, compared with that in the busbar boundary conditions [|J. 
Of much interest here are the recent results on the dynamical exponent. Under the pe- 
riodic boundary conditions, remarkably different values of the dynamical exponent have 
been obtained in the resistively shunted Josephson junction (RSJ) dynamics and in the 
time-dependent Ginzburg-Landau (TDGL) dynamics HH, suggesting that the two kinds 
of dynamics may be inherently different. In particular the dynamical exponent in the RSJ 
dynamics turns out to be much smaller than expected and apparently inapplicable to the 
flux-noise experiment |7| . On the other hand, no difference has been observed when the fluc- 
tuating twist boundary conditions (FTBC), forcing the vortex interaction to be periodic ||, 
is employed ||; this implies that the difference in the dynamical exponent stems from the 
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employed boundary conditions rather than from the intrinsic property of the dynamics. 

This paper is to clarify the boundary effects and thus to resolve the discrepancy between 
the two types of the boundary conditions. For this purpose, we investigate in detail the 
current- volt age (IV) characteristics and the corresponding dynamic critical behaviors of 
the system in the CBC and in the FTBC. The finite-size-scaling analysis is performed to 
reveal that the two types of the boundary conditions yield different values of the dynamical 
exponent. In particular, the data obtained in the CBC are found not to display proper 
scaling behavior. We show that the discrepancy has its origin in the rather large boundary 
effects in the CBC and give a simple proposal that the proper scaling analysis in the CBC 
should be performed with the boundary data discarded. It turns out that the resulting IV 
characteristics indeed exhibit good scaling behavior and are consistent with those from the 
FTBC, confirming that the poor scaling behavior originates from the peculiar boundary 
effects in the CBC and that correct description of the dynamic behavior can be achieved by 
appropriate extraction of the bulk properties from the data. 

We begin with the set of equations of motion describing the dynamics of an RSJ array 
with critical current I c and shunt resistance R s : 

= ir- (i) 

Here faj = (pi — (pj is the phase difference across a junction with fa being the phase of the 
superconducting order parameter on grain i (at position r$), If xt is the external current fed 
into grain i, and the primed summation runs over the nearest neighbors of grain i. The 
thermal noise current Ly is assumed to be a white noise satisfying 

(L i:i (t + r)L kl (t)) = -|_<J(r)(5 ife fy - 6«6 jk ) (2) 

with the angular bracket denoting an ensemble average and T being the temperature. In 
the CBC, uniform current I is fed into each grain at one edge (x = 0) and extracted from 
each grain at the other edge (x = N), leading to Zf xt = I(S Xi o — S Xi x) |J. The boundary 
conditions for the phase variables are taken to be free in the x direction and periodic in the 
y direction. The average voltage across each junction is given by the relation 
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2eR q dt 



+ L sin fat + L, 



where (. . .) t stands for the time average. In the FTBC, on the other hand, periodic boundary 
conditions are imposed on {</>j} in both directions, yielding the equations of motion in the 
form || 

3 s 

with Tij = Ti — Tj. The dynamics of the twist variables A = (A^., A y ) is governed by 

H dA x I c I 

^sia^ij - A y ) + L Ay , (5) 
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where J2{ij) x , v s denotes the summation over all nearest-neighboring pairs in the x(y) direction 



and the thermal noise current La , , satisfies 



2k T 

(L^ ) {t + T)L± w {t)) = J £ r 6{T). (6) 



The average voltage can then be obtained from the time evolution of the twist variable A x : 

To obtain the IV characteristics in the two boundary conditions, we integrate directly 
the coupled equations of motion given by Eq. (JJ) for the CBC and Eqs. (|J) and @ for the 
FTBC, with the time step At = 0.05 (in units of H/2eR s I c ). The data are averaged over 30 
independent runs while in each run the averaging time is chosen in such a way that at least 
25 vortices are created. We concentrate on the temperature T = 0.84 (in units of hl c /2e), 
which is below the BKT transition temperature T c w 0.89 |I0]| , and henceforth write the 
current and the voltage in units of I c and I C R S , respectively. 

Figure [I] presents the resulting IV characteristics in the CBC and in the FTBC. We 
have considered several system sizes ranging from N = 4 to iV = 64, to clarify size effects. 
At high currents the average voltage drop is found to be almost independent of the size 



in both boundary conditions. As the current is decreased, however, the system undergoes 
a crossover from the nonlinear-resistance regime to the linear-resistance regime, exhibiting 
size-dependence. The crossover current apparently decreases with the increase of the size, 
suggesting that the linear resistance at low currents reflects the finite-size effects PJTTf rather 



than the lattice effects ||. Here one can observe two striking differences between the IV 
characteristics in the CBC and those in the FTBC: First, at low currents the voltage drop in 
the CBC is much higher than that in the FTBC. Furthermore, the ratio of the former to the 
latter grows larger as the system size is increased. Secondly, the decrease of the crossover 
current with the system size is faster in the FTBC than in the CBC; this implies that the 
FTBC is more effective in reducing the finite-size effects. 

We now analyze the IV data by means of the finite-size scaling method, and investigate 
the dynamic critical behavior both in the CBC and in the FTBC. We assume the dynamic 



scaling form Hi] 



NR 1/z = f(NI), (8) 

where R is the resistance of the system and z is the dynamical exponent. The dimensionless 
scaling function f(x) is expected to possess the asymptotic behavior: f(x) oc constant for 
small x and f(x) oc x for large x. The small-x behavior follows from the behavior of the 
linear resistance for / — ► 0: 

R L ~ N~ z (9) 

whereas the finiteness of the resistance in the limit N — > oo elucidates the behavior of f(x) 
for large x. 

We first fit the linear resistance Rl in the CBC to the form in Eq. (£5|), and display in 
the inset of Fig. 0(a) the log-log plot of the linear resistance versus the system size, which 
shows that the linear resistance is well described by Eq. (|9]) with the dynamical exponent 
z = 0.91 ± 0.07. Note, however, that this is presumably deceptive since the value of z 
less than two is likely to be erroneous, especially, below the transition temperature: In a 



two-dimensional Josephson-j unction array, the value is generally believed to be z = 2 at 



the transition temperature and to grow larger as the temperature is lowered pj| . Indeed 
the scaling plot with the value z = 0.91, which is presented in Fig. 0(a), manifests that the 
system in the CBC does not follow the one-parameter dynamic scaling given by Eq. (^). 

We next examine the IV data in the FTBC, and fit the corresponding linear resistance 
Rl- This yields the dynamical exponent z = 2.73 ± 0.02, as shown in the inset of Fig. @(b), 
which is quite different from the value in the CBC and indeed consistent with the results of 
previous studies [p,|Tl|. In particular, unlike the result of the CBC, the scaling plot shown in 



Fig. ^(b) reveals that the scaling form in Eq. @ is observed well in the FTBC. Namely, the 
IV data in the nonlinear-resistance regime as well as in the linear-resistance regime collapse 
nicely into a single curve and the scaling function follows the expected asymptotic behaviors. 

Both the values of the dynamical exponent and the resulting scaling plots indicate that 
the CBC are not appropriate for describing the dynamical behavior of the system whereas 
the FTBC give a correct description. This discrepancy is likely to have its origin in the effects 
of the boundary, which is formed for the purpose of the current injection in the CBC but is 
absent in the FTBC. To find a clue to the boundary effects, we examine the distribution of 
the voltage drop along the x direction, which is, according to Eq. (|3|), proportional to the 
gradient of the phase profile. The phase profile 

M*) = i EIW*-) - <Wo)L (io) 

y 

where r is the average time, has been obtained for several system sizes and currents, which 
turns out to be essentially the same as the profile obtained in Ref. |4[]. It demonstrates 
that the voltage drop across a few columns near the boundary is comparable to, or even 
larger than, that across the rest of the system. Such tendency becomes prominent at lower 
currents, and leads to an overestimation of the voltage drop in the system, especially at 
low currents; this may have significant effects on the dynamic behavior of a system of finite 
size |§. 

To make clear that the boundary effects are the main origin of the failure of the dynamic 
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finite-size scaling in the CBC, we consider the (iV+10) x N array, consisting of iV+10 columns 
and N rows under the CBC, and calculate the voltage drop between the grains at x = 5 and 
x = N+5. Namely, we propose that the proper scaling in the CBC should be performed 
with the boundary data discarded, which helps to estimate correctly the voltage drop in 
the bulk [T2|]. Here the number five has been chosen since the phase profile reveals that the 
boundary effects are severe, i.e., the anomalous voltage drop is substantial, across the five 
columns at each edge, apparently regardless of the system size. The IV characteristics for 
several sizes, obtained in this scheme, are plotted in Fig. || which shows marked difference 
from Fig. 0(a): The voltage drop is significantly reduced at small currents, resulting in the 
IV characteristics which look almost the same as those obtained from the FTBC. 

Indeed the finite-size scaling analysis, shown in the inset of Fig. f|, leads to the dynamical 
exponent z = 2.84 ± 0.06, which is consistent with that obtained from the FTBC. Further, 
with this value of the dynamical exponent, the TV data in Fig. |3] collapse remarkably well into 
a single curve, as illustrated in Fig. |j. ^From these results we thus reach the conclusion that 
the CBC give a correct description of the dynamical behavior of the Josephson-j unction array 
only when the rather large voltage drop near the boundary is discarded in the analysis. Such 
boundary effects should be irrelevant in the thermodynamic limit, but can affect significantly 
the properties of a finite system, generating substantial errors in the scaling analysis. Taking 
into account that the periodic boundary conditions are naturally generalized to the CBC 
in the presence of driving currents, we can infer that the erroneous value of z obtained in 
Ref. || also has its origin in the boundary conditions rather than in the RSJ dynamics itself. 

In summary, we have studied the boundary effects on the dynamical behavior of a two- 
dimensional RSJ array, both in the CBC and in the FTBC. The finite-size-scaling analysis 
has been performed to reveal not only the discrepancy in the value of the dynamical exponent 
but also the lack of scaling in the CBC, steming from the rather large boundary effects. We 
have thus given a simple proposal that in the CBC the proper scaling should be performed 
with the boundary data discarded, and shown that such removal of the boundary effects 
indeed restores the correct scaling behavior. This confirms the crucial role of boundary 
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conditions in the finite-size scaling analysis and resolves the discrepancy between the two 
types of the boundary conditions. 
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FIG. 2. Scaling plots of the IV data (a) in the CBC with z = 0.91 and (b) in the FTBC 
with z = 2.73. Each inset displays the log- log plot of the linear resistance versus the size and its 
power-law fitting curve, from which the value of the dynamical exponent z is estimated. 
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FIG. 3. IV characteristics of (N+10)xN square Josephson-j unction arrays at T = 0.84 for 
several sizes in the CBC. The voltage is measured between the grains at x = 5 and at x = iV+5. 
Lines are merely guides to the eye. 
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FIG. 4. Scaling plot of the IV data presented in Fig. ||, with z = 2.84. Inset: the log-log plot 
of the linear resistance versus the size and its power-law fitting curve. 
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